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ABSTRACT 

We present a survey for bright Kuiper Belt Objects (KBOs) and Centaurs, 
conducted at the Kitt Peak National Observatory (KPNO) 0.9 m telescope with 
the KPNO 8k Mosaic CCD. The survey imaged 164 sq deg near opposition to 
a limiting red magnitude of 21.1. Three bright KBOs and one Centaur were 
found, the brightest KBO having red magnitude 19.7, about 700 km in diameter 
assuming a dark Centaur-like 4% albedo. We estimate the power-law differential 
size distribution of the Classical KBOs to have index q = 4.2lo;3, with the 
total number of Classical KBOs with diameters larger than 100 km equal to 
4:.7t\i X 10^ Additionally, we find that if there is a maximum object size in 
the Kuiper Belt, it must be larger than 1000 km in diameter. By extending our 
model to larger size bodies, we estimate that 30^12 Charon-sized and 'i.2~^ij 
Pluto-sized Classical KBOs remain undiscovered. 

Subject headings: comets: general — Kuiper Belt, Oort Cloud — solar system: 
formation 
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1. Introduction 

The region beyond Neptune is populated by ~ 10^ Kuiper Belt Objects (KBOs) 
with diameter D > 100 km, with a total mass of ~ 0.2 Earth masses (Jewitt, Luu, & 
Trujillo 1998). These bodies are found in three dynamical classes: (1) the Classical KBOs, 
which have semimajor axes in the 41 AU < a < 47 AU range with low eccentricities, 
e < 0.15; (2) the Resonant KBOs, which occupy the mean-motion resonances with Neptune, 
predominantly the Plutinos (a ^ 39.4 AU) and the 2:1 objects (a ~ 47.7 AU); and (3) the 
Scattered KBOs, which are distant (a > 50 AU) and have highly eccentric orbits (e ~ 0.6; 
Jewitt & Luu 2000). To date, very little is known about the physical properties of these 
numerous Kuiper Belt Objects (KBOs) because most of the known objects are quite faint; 
the known KBOs have median red magnitude mn ^ 23.0. Studies of other outer Solar 
System objects suggest that bright {rriR < 20) bodies are ideal targets for physical studies 
(Brown et al. 1997; Luu & Jewitt 1998; Brown, Cruikshank & Pendleton 1999; Luu, Jewitt 
& Trujillo 2000; Brown 2000). They may also allow (1) albedo measurements through 
the use of combined thermal and visible measurements, as collected from Centaur 10199 
Chariklo (provisional designation 1997 CU26; Jewitt & Kalas 1998; Altenhoff, Menten & 
Bertoldi 2001) and KBO 2000 WRioe (Jewitt & Aussel 2001); and (2) direct imaging with 
the Hubble Space Telescope, such as performed on Pluto (Stern, Buie & Trafton 1997) and 
Centaur 2060 Chiron (Parker et al. 1997). Large KBOs offer higher chances for occultations 
and longer duration events, such as have been observed for Chiron (Elliot et al. 1995; Bus 
et al. 1996), Charon (Walker 1980; Elhot & Young 1991), and Pluto (Elliot et al. 1989; 
Milhs et al. 1993). 

Very few surveys have been sensitive to these bright KBOs, primarily due to the fact 
that the surface density of ecliptic KBOs at tur ~ 20 is low, about 1 per 100 sq deg (Ferrin 
et al. 2001; Larsen et al. 2001; and Trujillo, Jewitt & Luu, 2001). In the past, photographic 
plate surveys have been used to search for these objects. Tombaugh (1961) examined about 
20,000 sq deg in a decade-long survey to visual magnitude my ~ 15.5, and found only Pluto 
in 1930. Kowal (1989) searched 6400 sq deg to visual magnitude my ~ 20 and found one 
Centaur, 2060 Chiron, but no KBOs. Although plates easily image large areas, they suffer 
from poor sensitivity and surface defects. In addition, computer-assisted analysis of plates 
is more difficult than that of charge-coupled devices (CCDs). More recently, CCDs, which 
have a much higher quantum efficiency than photographic plates, have played a growing 
role, due to technological advancements allowing increased detector size. For a review of all 
published wide- field (> 50 sq deg) CCD surveys, see Table |l|. 

In this work we report results from a new survey of the outer solar system undertaken 
at the Kitt Peak National Observatory (KPNO) using a large format Mosaic CCD array 
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(Muller et al. 1998). We have surveyed 164 square degrees to limiting red magnitude 
rrifi^Q = 21.1 (defined as the magnitude at which survey efficiency drops to 1/2 the 
maximum value). This survey yielded 3 bright KBOs and 1 Centaur, the brightest KBO 
being mn = 19.7, corresponding to a diameter of 710 km, assuming a 4% red albedo. We 
use these data to constrain the number of large KBOs and place limits on the maximum 
size of the KBOs using a realistic model of the discovery process. 

2. Survey Data 

Observations were made with the now decommissioned KPNO 0.9 m f/7.5 telescope 
and the 8192 x 8192 pixel Mosaic CCD array. This array is composed of eight 4096 x 2048 
pixel CCDs, each with an independent bias level and flat-field response. Observations were 
made through a standard Johnson R ffiter to enhance the detection of the KBOs, which 
have V — R ^ 0.6. The plate scale was 0.43 "/pixel for the 15 /xm pixels, corresponding to 
a 0.96 sq deg field of view. This telescope and camera combination was unique in providing 
a nearly 1 sq deg field of view while simultaneously Nyquist-sampling stellar images. In 
addition, the chips were relatively free from defects and had typical quantum efficiencies of 
~ 0.85 in R. Telescope and detector parameters are listed in Table ^. Fields were chosen to 
be within 1.5 hours of opposition where the parallactic motion of the KBOs is greatest (~ 3 
"/hr), and easily distinguishable from the main-belt asteroid motion (;:^ 25 " /hi). Nearly 
all (~ 90%) of the fields were confined to be within 5 deg of the ecliptic, as depicted in 
Figures ^ and |^. Each field was imaged three times, using 200 s exposures with a ~ 1 hour 
separation between images. The observed fields included in the survey passed two quality 
tests: (1) they must have been taken during photometric conditions, and (2) they must 
have had a characteristic stellar Full- Width at Half-Maximum (FWHM) better than 2.5 ". 
Of the 285 fields imaged in the survey, only 171 passed these criteria. These fields are listed 
in Table |^. As detailed in Table |^, 1998 SNigs was discovered by Spacewatch three nights 
before we serendipitously detected the object. We verified that no known objects remained 
undiscovered in our data by computing the ephemerides of all bright (my < 22.5) objects 
with opposition motion slower than 10 arcsec/hr for the epochs observed. 

The width, ellipticity, and position angle of the Point Spread Function (PSF) varied 
significantly over the Mosaic CCD field of view. We quantified this variation by fitting an 
elliptical Moffat PSF (Moffat 1969) to stars in the linear flux regime of the CCDs for a large 
number of images. Figure ^ shows the variations in the FWHM of the PSF minor axis, 
as measured by IRAF's IMEXAM procedure. The ellipticity was found to be correlated 
with the minor axis FWHM, ranging in magnitude from 0.1 - 0.3 over the focal plane. 
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The width, ellipticity and position angle information was used to add simulated KBOs to 
the data in order to test our Moving Object Detection Software (MODS, Trujillo & Jewitt 
1998). 

The seeing in the survey varied from 0.9"to 2.4", with median 1.5". We tested MODS 
with 3 groups of images based on the seeing (0.9"- 1.4", 1.4"- 1.9", and 1.9"- 2.4"; with 
median seeing 1.3", 1.6", and 2.1 ", respectively). Artificial moving objects, with profiles 
matched to the characteristic PSF for each image group and location in the focal plane, 
were added to flattened images. A tally of objects was recorded, and the detection efficiency 
(Figure ^ was computed and found to be uniform with respect to sky-plane speed in the 
1.5 - 10 "/hr range. 

The PSF- and seeing-corrected efficiency function was fitted with the same hyperbolic 
tangent function used in Trujillo, Jewitt & Luu (2000), given by 

. = ^ (tanh ( "^^Y"^^ ) + l) . (1) 

The limiting red magnitude ruji^Q is the brightness where e = emax/2, with Smax defined as 
the maximum detection efficiency obtained for bright objects. The characteristic magnitude 
range over which the efficiency drops from ^max to zero is a. The values of these parameters 
for each seeing range are listed in Table ^. The poorer seeing images show a drop in m^^so 
compared to those with good seeing, consistent with reduced signal-to-noise ratio. Thus, 
the limiting magnitude can be characterized as m^so = 21.8 — 2.5 log(^), where 6 represents 
the FWHM of a stellar image in arc sec. The quantity ^max varies little for the different 
seeing categories; even for the worst seeing, ^max is reduced by only 9% from the best seeing. 

The discovery circumstances for the 4 KPNO objects appear in Table ^. Although the 
discovered objects were found during good seeing (Table H), the null detection of objects in 
the medium and poor seeing cases is not statistically significant. Given the a = 0.66 slope 
of the Cumulative Luminosity Function (§ ^) and the limiting magnitude and sky area 
presented in Table |, 46% of the objects should be found during good seeing. Thus, the 
probability of all 4 object detections occurring by chance during good seeing is 0.46^ = 4.4% 
(less probable than a 2a Gaussian event). We therefore combine the data for all seeing 
cases into a single global efficiency function. 

3. The Kuiper Belt 

The four critical quantities we estimate for the KBOs are: (1) the Cumulative 
Luminosity Function, or CLF, (2) the size distribution and total number of objects, (3) the 
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KBO maximum size, and (4) the total number and mass of large KBOs 2000 km and 
^ 1000 km diameter). 

3.1. The Cumulative Luminosity Function 

The Cumulative Luminosity Function (CLF) describes the surface density of KBOs 
near the ecliptic (S) versus limiting red survey magnitude {niR) and is fit by the equation 

log = a{mR- mo), (2) 

where a describes the slope and mo is the red magnitude at which E = 1 KBO per sq deg. 
Assuming heliocentric and geocentric distance and albedo independent of object size, and a 
differential size distribution n{r)dr oc r~'^dr, where n{r)dr describes the number of bodies 
with radii between r and r + dr, the size distribution of the KBOs can be directly measured 
via (Irwin et al. 1995) 

q = 5a + l. (3) 

We estimated the CLF by combining the discovery results of the 86 KBOs found by 
Trujillo, Jewitt & Luu (2001, hereafter OITJL) with the three KBOs found in our data. 
Although the OITJL data dominate the total numbers of faint objects, the two surveys 
have found comparable numbers of bright objects. In addition, our brightest object (1999 
DEg, mji = 19.7) was over a magnitude brighter than the brightest object from OITJL 
(1999 CDi58, mji = 21.0). These two data sets improve upon previous survey results due to 
discovery statistics that cover nearly a 5 magnitude range (19.7 < m/j < 24.4). 

Since the CLF is a cumulative measure, the CLF value at a given magnitude is 
correlated with the CLF values at fainter magnitudes. Thus, a linear fit to the CLF will 
implicitly weight the bright bodies (which are counted in the bright and faint data points) 
to a greater degree than the faint bodies (which are only counted in the faint data points). 
We instead perform a maximum-likelihood fit to find a and mo (Gladman et al. 1998), 
using the efficiency function (Equation [^) and the sky area imaged (164 sq deg). It was 
impossible for all fields to be centered on the ecliptic due to the large sky area imaged. 
Thus, sky coverage was computed by weighting each field imaged by the fraction of KBOs 
expected at that echptic latitude, assuming the best-fit inclination distribution of OITJL, 
a half-width of 20 deg. For our KPNO data, this correction reduces our 164 sq deg to an 
effective area of 95.7 sq deg, resulting in a minor change (23%) in the normalization of 
log(S), a factor 2 smaller than the la Poisson uncertainty for log(S) due to the 4 KBOs 
discovered (45%). The resulting best fit follows a = 0.66 ±0.06 and mo = 23.32 ±0.09. This 
result is in statistical agreement with many other previous works (Jewitt, Luu & Trujillo 
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1998 with a = 0.58 ± 0.05, mo = 23.27 ± 0.11; Gladman et al. 1998 with a = 0.76;^;}? and 
mo = 23.40tJ]:i^; Chiang & Brown 1999 with a = 0.52 ± 0.05 and mo = 23.5 ± 0.06; and 
OITJL a = 0.63 ± 0.06, mo = 23.04lo;o9). We present this simple fit, our observations, and 
other works covering > 50 sq deg in Figure ^ Using equation |^, we find the exponent of 
the differential size distribution to be g = 4.3 ± 0.3. We present a more detailed analysis in 
the next section, where we use a more realistic simulation to constrain the slope of the size 
distribution. 



3.2. The Size Distribution of the KBOs 

To better constrain the size distribution, we estimate q and the total number 
of Classical KBOs larger than 100 km, iVcKBOs(-D > 100 km), using a more realistic 
maximum-likelihood simulation than the CLF fit described in the previous section. Bias 
factors such as the detection efficiency and the distribution of discovery distances are 
correctly handled in this more detailed simulation. A full description of our simulation 
can be found in Trujillo, Jewitt & Luu (2000) and OITJL. We summarize the assumed 
quantities in Table ^ and key procedures here: 

(1) Simulated orbital elements and object sizes were drawn for the Classical KBOs 
(CKBOs). 

(2) Object brightnesses were computed from the following formula (Jewitt & Luu 
1995): 

m = m© - 2.5 log(p$(a')^^) + 2.5 log(2.25 x lO^^^i^^A^), (4) 

where a' represents the phase angle of the observations, $(q;') is the Bowell et al. (1989) 
phase function for dark bodies {G = 0.15), geometric albedo is given by p, object radius 
is described by r [km], and R [AU] and A [AU] describe the heliocentric and geocentric 
distance of the body, respectively. The apparent red magnitude of the Sun is ttiq = —27.1. 

(3) Object ecliptic coordinates and velocities were computed from the equations of 
Sykes & Moynihan (1996, sign error corrected). 

(4) The observed field areas and efficiency functions were used to determine which of 
the simulated objects could be detected in either our survey or that of OITJL, based on the 
ecliptic coordinates, velocities and brightnesses computed in steps 2 and 3. 

(5) Two tallies of the radii of the "detected" simulated objects were kept, one for our 
KPNO 8k Mosaic survey and one for the OITJL survey. Steps 1-4 were repeated until 
there were a factor 10 to 100 more simulated objects in each radius bin than were actually 
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observed. 

(6) The likelihood of the model producing the observed distribution of radii was 
estimated by assuming Poisson detection statistics with the number of "detected" simulated 
objects providing the expectation value for each bin. This dataset consists of the 57 CKBOs 
found by OITJL, and the single CKBO found in the KPNO survey (1999 DFg). Although 
the OITJL data contains many more objects, the KPNO survey has comparable area 
coverage, and number of bright object discoveries. 

(7) Steps 1-6 were repeated, each time varying the size distribution index, q, and the 
total number of Classical KBOs, A^^ckbosI-D > 100 km), with the ultimate goal of finding 
the number of CKBOs and size distribution that has the maximum likelihood of producing 
the observations. 

The primary process we are modeling is the effect of heliocentric distance R and 
geocentric distance A on the apparent red magnitude m^j, and the resulting probability 
of detection based on detection efficiency (Equation |1|) and survey area. The inclination 
distribution assumed in Table ^ is equivalent to the best-fit value found by OITJL — 
a Gaussian model with half-width of 20 deg. The true ecliptic latitude distribution of 
observed fields was used in the simulation, so no explicit correction for finite belt thickness 
is necessary, as needed in the CLF fit (§ lOI ). 

The results for the best-fit size distribution index, q, and number of Classical KBOs 
with diameters D > 100 km, A'ckbos(-D > 100 km), are presented in Figure |. This figure 
depicts the contours of constant probability that the given q and A'ckbos(-D > 100 km), 
could produce the observed distribution. We find that 

q = 4.2i[]-^ (la) and 

= 4.21J:° (3a), 
and 

A'cKBOs(£' > 100 km) = 4.7l|-[; x 10^ (la) and 

= A.rltix 10^ (3a). 

This estimate of q is in formal agreement with that derived from the CLF alone at < Icr 
level. 

The results for q to not depend sensitively on our model assumptions. Our values 
for A'cKBOs(-D > 100 km) are somewhat dependent on the inchnation distribution 
assumed, as thicker inclination distributions can hide more objects at high ecliptic 
latitudes where no observations were made. For example, increasing the half- width of the 
inclination distribution by a factor of two, from 20 deg to 40 deg results in an increase 
in the best-fit A'ckbos(^ > 100 km), from 4.7 x 10^ to 7.5 x 10^, a factor 1.6. The 
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best-fit value for A/'ckbosI-D > 100 km) is roughly linear with half-width ii/2, following 
A^CKBOs(^ > 100 km) ~ (1.9 + O.Uii/2) x 10^. 

3.3. The Maximum Size of Kuiper Belt Objects 

Adopting the best-fit q = 4.2, we simulated the CLF that would be found by our 
survey given several different maximum sizes for the Classical KBOs, with results appearing 
in Figure 0. The expected CLF from the Classical KBOs was multiplied by the observed 
KBO : Classical KBO ratio (89:58) to construct this graph, which makes the implicit 
assumption that the CLF of the Classical KBOs matches that of the KBOs as a whole. 
This assumption is reasonable, as a maximum-likelihood fit to the CLF for the Classical 
KBOs alone (a = 0.7lt°om) 

is statistically consistent with the fit to all observed populations 

(a = 0.66 ±0.06, §|3l|). 

Our analysis rules out the -Dmax = 500 km and -Dmax = 250 km diameter maximum 
sizes, as is expected since we found one body with diameter D ~ 700 km and another with 
D ~ 500 km. The Z^max = 1000 km simulation is the best-fit to our data, however, these 
data are formally consistent will all greater -Dmax- Thus, if there is a maximum size cutoff, it 
is greater than or equal to -Dmax = 1000 km, because we do not observe a significant surface 
density decrease for bright objects beyond that described by the linear slope of the CLF. 
This analysis is in agreement with all other published wide-field surveys except for Kowal 
(1989). Kowal's (1989) datum could be explained by a precipitous drop in surface density 
between 19 < < 20, as would be expected if there were a maximum size of -D = 1000 
km to the KBOs. Such a model would be in violation of the Tombaugh (1961) datum point, 
unless Pluto is considered a special case. In either case, our observations indicate that if 
there is a maximum size, it must be -Dmax > 1000 km. 

3.4. The Total Number and Mass of Large Kuiper Belt Objects 

The number of large KBOs is of fundamental interest because it is a direct measurement 
of the runaway accretion in the planet forming era of the solar system (Kenyon & Luu 
1998). Two sizes of typical interest are the number of Charon-sized bodies (~ 1000 km in 
diameter) and the number of Pluto-sized bodies (~ 2000 km in diameter). The number of 
bodies with diameters larger than -Dmm, given by N{D > -Dmin), can be directly calculated 
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from the results of our maximum likelihood simulation. For g 7^ 1, 



N{D > D^^) = iVcKBO.(^ > 100 km) 



100 km 



min 




(5) 



Thus, from Figure we find that the total numbers of Charon- 
in the Classical Kuiper Belt follow 



and Pluto-sized bodies 



N{D > 1000 km) = SOtH (Ict) and 
= 30tll (3a), 



and 



N{D > 2000 km) = 3.2lff (la) and 



Thus, assuming that the q = 4.2 size distribution continues to large sizes, we expect that 
there are > 1 {3a lower-limit) undiscovered Pluto-sized bodies in the Classical Kuiper Belt 
and > 8 {3a lower-limit) Charon-sized bodies, of which only a few have been discovered. 
If the size distributions of the Classical and Scattered KBOs are identical, then a similar 
number of bodies should be present in the Scattered Kuiper Belt, which has approximately 
the same number of bodies (Trujillo, Jewitt & Luu 2000). These results are consistent with 
the Kuiper Belt growth models of Kenyon and Luu (1998, 1999) using velocity evolution and 
coUisional fragmentation. In these models, several Pluto-sized bodies grow concurrently in a 
low-mass (~ IOMq) disk. These models have two characteristic observable parameters: (1) 
The size distribution of the resulting Kuiper Belt corresponds to g ~ 4, and (2) several large 
bodies are formed concurrently. Their favored model produces ~ 150 Charon-sized bodies 
and 1 Pluto-sized body in 36.5 Myr (Kenyon and Luu 1999a) and is roughly consistent with 
our observations, although no prediction is made about the subsequent destruction rates of 
such bodies from catastrophic collisions. 

The possibility of finding these Pluto- and Charon- sized bodies can be estimated 
by examining the discovery conditions of other bright bodies. Table ^ lists the three 
KBOs to date brighter than red magnitude rriR = 20 that have been found in published 
surveys (with the exception of Pluto). These three bodies were found in the heliocentric 
distance range 30-39 AU. Assuming opposition observations and 4% albedo, 2000 km and 
1000 km diameter bodies at these heliocentric distances would have have brightnesses 
16.9 < niR < 18.2 and 18.4 < < 19.7, respectively, according to Equation ^. These 
magnitudes are within range for most time-resolved wide-field survey telescopes, however 
the apparent motion of these bodies, ~ 3 "/hr, is considerably slower than ^ 25 "/hr 
main-belt asteroids. Since most such observing programs are designed to find near-Earth 
asteroids and other high proper motion planetary bodies, plate scales and timebases 



-10- 



between images are not usually conducive to finding the KBOs. Thus, special observing 

procedures or software methods must be adopted to find the slower KBOs, as has been 
done by Spaccwatch (Larsen et al. 2001). A time-resolved all-sky rriji = 20.0 magnitude 
survey should be able to find all the remaining Pluto- and Charon- sized bodies, as long as 
suitable measures are taken to detect ~ 3 "/hr bodies. Such a survey could be conducted 
with a dedicated 1 meter class telescope outfitted with a CCD camera allowing a wide-field 
of view, similar to the KPNO 0.9 m and 8k CCD Mosaic instrument combination used in 
this work. 

These analyses presume that the KBOs have a 4% albedo, consistent with the Centaurs 
(Jewitt & Luu 2000). If the mean KBO albedo, p, departs from this value, then the radial 
sizes of the bodies will be affected by a factor (p/0.04)~'^/^. For a 14% albedo such as that 
of Chiron (Campins et al. 1995) and a 38% albedo such as found for Charon (Buie, Tholen 
& Wasserman 1997), this amounts to a size factor 1.9 and 3.1, respectively. 

The total amount of mass in 1000 km to 2000 km diameter Classical KBOs for our 
best-fit q = 4.2 case is 



with p representing the bulk density, p denoting the albedo, and — 6.0 x 10^ kg. This 
accounts for 18% of the total mass found in 100 km < D < 2000 km bodies. Again, the 
Scattered KBOs may contain a similar amount of material. 



Although there is no formal International Astronomical Union definition of the Centaur 
population, we classify bodies whose perihelia fall between that of Jupiter and Neptune as 
Centaurs, similar to the Larsen et al. (2001) definition. One such Centaur was found in the 
KPNO 8k data (2000 EE173). We therefore estimate the sky-plane density of Centaurs to 
be l.OlJI X 10-2 (la) Centaurs per sq dcg at our survey limiting magnitude m^so = 21.1. 
This is consistent with other recent estimates of the Centaur population: Larsen et al. 
(2001) predict a value of 1 x 10"^ at our limiting magnitude and Sheppard et al. (2000) 
predict a similar surface density of 3 x 10^^ Centaurs per sq deg. 



M(1000 km < D < 2000 km) = 5.4 x 1O~^M0 



1000 kg m 



P 




4. 



The Centaurs 
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5. Summciry 

We have conducted a survey for KBOs and Centaurs with the KPNO 0.9 m telescope 
and 8k Mosaic CCD, which together provide a nearly 1 sq deg field of view, reaching red 
magnitude rriR = 21.1 in 200 s in typical seeing conditions at Kitt Peak. We detected 3 
KBOs and 1 Centaur in 164 sq deg examined near the ecliptic. Combining these new data 
with the results of Trujillo, Jewitt & Luu (2001), we find the following results: 

(1) The Cumulative Luminosity Function slope is a = 0.66 ± 0.06, with 1 KBO / sq 
deg occurring at red magnitude mo = 23.32 ± 0.09. 

(2) The Classical KBOs follow a g = 4.21q;3 (la) differential size distribution, 
with the total number of Classical KBOs with diameters larger than 100 km given by 
-^CKBOs(-D > 100 km) = 4.7lJ;o ^ 10"^) assuming a Gaussian inclination distribution with 
half-width 20 deg. 

(3) No statistically significant evidence for a maximum size cutoff to the KBOs was 
found. However, if such a maximum diameter exists, it must have i^max > 1000 km. 

(4) Assuming the differential size distribution extends to large sizes, the total number 
of Classical KBOs with diameters D ;^ 1000 km and D ;^ 2000 km is 30t\t and 3.2ii:?, 
respectively (Icr). The bodies in the 1000 km < D < 2000 km diameter range account 
for ~ 1/5 of the total mass in the 100 km < D < 2000 km diameter range, representing 
M(1000 km < L> < 2000 km) = 5.4 x 1O-^M0, assuming a bulk density of 1000 kg m'^. 

(5) The sky-plane surface density of Centaurs brighter than ttir — 21.1 is I.OIqIs ^ 10~^ 
(Icr) Centaurs/sq deg near the ecliptic. 

We thank the National Optical Astronomy Observatories (NOAO) for granting 
telescope time, David Tholen for computation of the orbital elements, Amanda Sickafoose 
and Wyn Evans for observational assistance, David Jewitt for editorial comments, and 
Nigel Sharp for technical assistance at the 0.9m telescope. This work was supported in part 
by a NASA grant to David C. Jewitt, NASA grant NAG5-3940 to JLE and through NSF 
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Fig. 1. — Fields that passed quality tests, imaged during Sep 1998. The ecliptic is denoted 
by a solid line. 
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Fig. 2. — Fields that passed quality tests, imaged during Feb 1999 (white squares) and Mar 
2000 (grey squares). The ecliptic is denoted by a solid line. 
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Fig. 3. — Minor axis Moffat FWHM in pixels (1 pixel = 0.43 ") of a stellar profile across 
the Mosaic CCD. Note that the chip is misaligned with the optical axis of the telescope by 
about 1500 pixels = 2.3 cm. 
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Fig. 4. — The efficiency function of the 8k CCD on the KPNO 0.9m in 1.3"seeing, with a fit 
to the hyperbohc tangent efficiency function. 
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Red Mognitude 



Fig. 5. — Our measurement of the Cumulative Luminosity Function (CLF), which represents 
the number of KBOs deg"^ near the echptic (filled circles) brighter than a given apparent red 
magnitude. The line represents the best fit {a = 0.66 ± 0.06 and mo = 23.32 ± 0.09) to the 

CLF determined by this survey and that conducted by Trujillo, Luu & Jcwitt (2001, OITJL). 
Other points arc from previous works: OIF ct al. (Fcrrin et al. 2001); OIL et al. (Larsen 
et al. 2001); OOS ct al. (Sheppard et al. 2000); 98JLT (Jewitt, Luu & Trujillo 1998); 89K 
(Kowal 1989); 88LJ (Luu & Jewitt 1988); and 61T (Tombaugh 1961), with arrows denoting 
upper limits. 
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Fig. 6. — The normalized probability that the observed distribution of radii match the 
expected distribution given the free parameters: (1) g, the size distribution index, and (2) 
A^CKBOs(-D > 100 km), the total number of Classical KBOs (sigma values correspond to 
Gaussian confidence limits: la — 68.3%, 2a — 95.4%, etc.). The maximum hkelihood occurs 
at? = 4.2l°:|(l(7). 
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Fig. 7. — The expected observed CLF for the q = 4.2 case, given maximum KBO diameters 
(left-to-right) of 12,000 km (~ Earth-sized), 1000 km (~ Pluto-sized), 500 km Charon- 
sized) and 250 km. The 1000 km case represents the best fit, however only the 250 km and 
500 km cases are rejected at the 3a level. Formally, there is no evidence for an upper-size 
cutoff. 
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Table 1. Wide-Field CCD Surveys sensitive to KBOs 



Area 




Centaurs'' 


KBOs^ 


Brightest KBO 


Source 


[sq deg] 
















1428 


18.8 R 


0(1) 


0(0) 








Sheppard et al. (2000) 


522 


20.5 V 


3(1) 


7(5) 


1999 GQ21 


[mv = 


20.5) 


Larsen et al. (2001) 


164 


21.1 R 


1 (0) 


2(1) 


1999 DEg 


{mR = 


19.7) 


This Work 


73 


23.7 R 


0(0) 


84 (2) 


1999 CDi58 


{rriR = 


21.0) 


Trujillo, Jewitt & Luu (2001) 


67 


20.1 R 


0(0) 


1(0) 


2000 EB173 


ijUR = 


19.3) 


Ferrin et al. (2001) 


52 


22.5 R 


0(0) 


13 (0) 


1996 TL66 


{ruR = 


20.6) 


Jewitt, Luu & Trujillo (1998) 



Note. — Published CCD surveys with > 50 sq deg have been included. 
^Limiting magnitude of survey. 

^The number of Centaurs discovered is listed, with serendipitous detections of previously discovered 
objects listed in parentheses. 

•^The number of KBOs discovered is listed, with serendipitous detections of previously discovered 
objects listed in parentheses. 
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Table 2. KPNO Survey Parameters 



Quantity KPNO 0.9m 

Focal Ratio f/7.5 
Instrument KPNO 8k x 8k Mosaic 

Plate Scale ["/pixel] 0.43 

Field Area [deg^] 0.957 

Total Fields 171 

Total Area [deg^] 164 

mRso'^ 21.1 

6lf' ] 0.9-2.4 

Filter R 

Quantum Efficiency 0.85 



^the red magnitude at which detection 
efficiency reaches half of the maximum efficiency 

''range of minor axis Full Width at Half 
Maximum of stellar sources near focal plane 
center 
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Table 3. CFHT Field Centers 



ID 


UT date 


UT times 




h 

a ° 


6 ^ 


^ Objects 




1998 


Sep 


25 


2:54 


3:58 4:54 


1 9 


99-1 


09-99-01 


P 




1998 


Sep 


25 


3:00 


4:04 5:00 


in 

-LU 


99-9n-^/t 




P 


prim HQ 


1998 


Sep 


25 


3:06 


4:10 5:07 


q 


99-9^n6 




P 




1998 


Sep 


25 


3:13 


4:16 5:13 


a 
o 


99-9'^-l R 

ZjZj .ZjO.±0 


U-L .^c^.U 1 


P 




1998 


Sep 


25 


3:19 


4:22 5:19 


u 


99.97.Q1 


0*^-1 9-4Q 

uo . XZ .'ity 


P 




1998 


Sep 


25 


3:25 


4:29 5:25 


4 


99-9Q-ZL4 


04-'^fi-9Q 


P 




1998 


Sep 


25 


3:31 


4:35 5:31 


3 


22-31 -56 


-06-00-00 

\J\J •\J\J •\J\J 


m 


prim 09 


1998 


Sep 


25 


3:38 


4:41 5:38 


1 


99-34-1 1 


-07-38-On 

U 1 .oo.uu 


111 


prim nn 


1998 


Sep 


25 


3:44 


4:48 5:44 


n 


99-3fi-9fi 

£j £j .\J\J . £j\J 


-08-47-On 

UO.t: t .UU 


111 


ppIO? do 


1998 


Sep 


25 


5:50 


6:47 7:50 


n 

u 


Z(O.UO.(JU 


n'^-9R-00 
-uo.zo.uu 


P 


cClU 1 .Ul 


1998 


Sep 


25 


5:56 


6:53 7:56 


9 


ZiO.UU.OO 


OzL-0^-00 
-u^.uo.uu 


m 


ct^iU / .UO 


1998 


Sep 


25 


6:03 


7:00 8:03 


Q 
O 


9Q.nzL-l 7 


09-49-00 

-UZ .4:Z . UU 


m 


cCiUO.UO 


1998 


Sep 


25 


6:09 


7:06 8:09 


Q 
O 


9Q-nQ-'^Q 
ZiO.uy .Oc/ 


09-OQ-OO 

-UZ.Uc/.UU 


P 




1998 


Sep 


25 


6:15 


7:12 8:16 


9 


9Q.1 1 .r;Q 

ZiO. 1 J. .Ot/ 


O^'^l -00 
-UO. 01 .UU 


m 


ppIoj^ no 


1998 


Sep 


25 


6:22 


7:25 8:22 


u 


9Q.14.9n 


-U'y:.u^.uu 


P 


pHoq do 


1998 


Sep 


25 


6:29 


7:31 8:28 





23- 19-43 


-04-20-On 




prinQ ni 


1998 


Sep 


25 


6:35 


7:37 8:35 


2 


23-1 7-99 


-02-57-00 




prlDQ 03 

C^iUc/ .U(J 


1998 


Sep 


25 


6:41 


7:44 8:41 


O 


93-1 5-ni 


-01 -34-on 

U X .Ot: . UU 


P 


<=>ni9^ no 


1998 


Sep 


25 


8:50 9:50 10:53 


n 

u 


UU.'4:U.Uc/ 




P 




1998 


Sep 


25 


8:56 9:56 10:59 


9 




OR-90-1 ^ 
uu.zu. lu 


P 


ppiocc no 


1998 


Sep 


25 


9:03 10:10 11:05 


Q 
O 


OO-zLl -9^ 


U / .fiO.UfJ: 


P 


ecl26.03 


1998 


Sep 


25 


9:10 10:17 11:11 


3 


00:46:56 


08:17:41 


P 


ecl26.01 


1998 


Sep 


25 


9:16 10:23 11:17 


1 


00:49:17 


06:54:48 


P 


ecl26.00 


1998 


Sep 


25 


9:22 10:29 11:23 





00:51:38 


05:31:52 


P 


ecl27.00 


1998 


Sep 


25 


9:28 10:35 11:29 





00:57:08 


06:06:12 


P 


ecl27.01 


1998 


Sep 


25 


9:37 10:41 11:35 


1 


00:54:47 


07:29:13 


P 


ecl27.03 


1998 


Sep 


25 


9:43 10:47 11:41 


3 


00:52:26 


08:52:12 


P 


ecl02.12 


1998 


Sep 26 


2:36 3:33 4:29 


12 


22:23:58 


02:53:28 


m 


ecl02.11 


1998 Sep 26 


2:42 3:39 4:37 


11 


22:26:12 


01:29:55 


m 


ecl02.09 


1998 Sep 26 


2:48 


3:45 4:43 


9 


22:28:25 


00:06:19 


g 


ecl02.08 


1998 Sep 26 


2:54 3:51 4:50 


8 


22:30:39 


-01:17:15 


g 


ecl02.06 


1998 Sep 26 


3:00 3:57 4:57 


6 


22:32:52 


-02:40:49 


m 
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Table 3 — Continued 



ID 


UT date 


UT times 




a ^ 


S ^ 




Objects 


ecl02.05 


1998 Sep 26 


3:06 4:04 5:03 


5 


22:35:07 


-04:04:21 


m 




ecl02.03 


1998 Sep 26 


3:13 4:10 5:09 


3 


22:37:22 


-05:27:55 


g 




ecl02.02 


1998 


Sep 26 


3:19 4:17 5:15 


1 


22:39:37 


-06:51:22 


m 




CC102.00 


1998 Sep 26 


3:26 4:23 5:35 





22:41:52 


-08:14:50 


m 




ecUO.OO 


1 998 




26 


5:42 6:38 7:43 





23:25:06 


-03:45:45 


m 




ecllO.Ol 


1998 


Sep 


26 


5:48 6:44 7:49 


2 


23:22:44 


-02:23:02 


m 




ecllO.03 


1998 


Sep 


26 


5:54 7:00 7:56 


3 


23:20:23 


-01:00:15 


g 




eclll.03 


1998 


Sep 


26 


6:00 7:06 8:00 


3 


23:25:44 


-00:25:49 


g 




ecUl.Ol 


1998 


Sep 


26 


6:07 7:12 8:08 


2 


23:28:06 


-01:48:33 


g 


1998 SNi65 


ecUl.OO 


1998 


Sep 


26 


6:13 7:18 8:14 





23:30:29 


-03:11:16 


g 




ecll2.00 


1998 


Sep 


26 


6:19 7:25 8:20 





23:35:52 


-02:36:35 


g 




ecll2.01 


1998 


Sep 


26 


6:25 7:31 8:27 


2 


23:33:29 


-01:13:53 


m 




ecll2.03 


1998 


Sen 


26 


6:31 7:37 8:33 


3 


23:31:07 


00:08:46 


m 




ecl28.00 


1 998 






8:39 9:49 10:50 





01:02:37 


06:40:22 


m 




ecl28.01 


1 QQR 


Sep 


zu 


8:45 9:56 10:56 


1 


01:00:18 


08:03:29 


g 




CC128.03 


1998 


Sep 


26 


8:59 10:02 11:02 


3 


00:57:58 


09:26:32 


m 




CC129.03 


1998 


Sep 


26 


9:05 10:08 11:08 


3 


01:03:31 


10:00:40 


m 




ecl29.01 


1998 


Sep 


26 


9:12 10:14 11:15 


1 


01:05:50 


08:37:32 


m 




ecl29.00 


1998 


Sep 


26 


9:18 10:20 11:22 





01:08:09 


07:14:19 


m 




ecl30.00 


1998 


Sep 


26 


9:24 10:26 11:28 





01:13:41 


07:48:03 


m 




ecl30.01 


1998 


Sep 


26 


9:31 10:32 11:34 


1 


01:11:24 


09:11:22 


m 




ecl30.03 


1998 


Sep 


26 


9:37 10:38 11:40 


3 


01:09:05 


10:34:40 


g 




ecl31.00 


1998 


Sep 


26 


9:44 10:44 11:46 





01:19:15 


08:21:34 


g 




ecl03.06 


1998 


Sep 


27 


2:26 3:26 4:36 


6 


22:38:15 


-02:08:29 


m 




ecl03.05 


1998 


Sep 


27 


2:32 3:33 4:42 


4 


22:40:30 


-03:31:54 


m 




ecl03.03 


1998 


Sep 


27 


2:38 3:39 4:48 


3 


22:42:44 


-04:55:20 


m 




ecl03.02 


1998 


Sep 


27 


2:44 3:51 4:56 


2 


22:45:01 


-06:18:43 


P 




ecl03.00 


1998 


Sep 


27 


2:50 3:57 5:02 





22:47:18 


-07:42:04 


m 




ecl04.00 


1998 


Sep 


27 


2:56 4:04 5:08 





22:52:44 


-07:08:59 


m 




ecl04.02 


1998 


Sep 


27 


3:02 4:10 5:14 


1 


22:50:26 


-05:45:46 


m 




ecl04.03 


1998 


Sep 


27 


3:08 4:16 5:20 


3 


22:48:09 


-04:22:30 


m 




ecl04.05 


1998 


Sep 


27 


3:14 4:22 5:26 


5 


22:45:52 


-02:59:10 


m 
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Table 3 — Continued 



ID 


UT date 


UT times 


(3 ^ 




5 " 


e Objects 


ecl04.06 


1998 Sep 27 


3:20 4:30 5:33 


6 


22:43:37 


-01:35:52 


P 


ecllS.OO 


1 QQ8 




97 


6:22 7:18 8:25 





23:41:14 


-02:01:50 


P 


ecll3.02 






97 


6:28 7:24 8:31 


1 


23:38:52 


-00:39:11 


P 


ecll3.03 






27 


6:36 7:30 8:37 


3 


23:36:29 


00:43:26 


P 


ecll4.03 


1998 


Sep 


27 


6:41 7:36 8:43 


3 


23:41:51 


01:18:16 


P 


ecll4.02 


1998 


Sep 


27 


6:47 7:43 8:49 


2 


23:44:14 


-00:04:20 


P 


ecll4.00 


1998 


Sep 


27 


6:53 7:49 8:55 





23:46:37 


-01:26:56 


P 


ecll5.00 


1998 


Sep 


27 


6:59 7:55 9:01 





23:51:59 


-00:51:59 


P 


ecll5.02 


1998 


Sop 


27 


7:05 8:02 9:07 


1 


23:49:37 


00:30:35 


P 


ecll5.03 


1998 


Sep 


27 


7:11 8:09 9:13 


3 


23:47:14 


01:53:11 


m 


ecl31.01 


1 998 




27 

Zj i 


9:28 10:23 11:17 


1 


01:16:59 


09:45:00 


P 


ecl31.03 


1998 


Sep 


27 


9:34 10:29 11:23 


3 


01:14:40 


11:08:23 


P 


ecl31.04 


1998 


Sep 


27 


9:42 10:35 11:36 


5 


01:12:21 


12:31:42 


P 


ecl31.06 


1 998 




27 


9:48 10:41 11:42 


6 


01:10:01 


13:54:55 


P 


ecl32.06 


1998 


Sep 


27 


9:53 10:47 12:13 


6 


01:15:39 


14:28:38 


P 


ecl32.05 


1998 


Sep 


27 


9:59 10:53 11:49 


5 


01:17:59 


13:05:18 


P 


ecl32.03 


1998 


Sep 


27 


10:05 10:59 11:55 


3 


01:20:17 


11:41:53 


P 


ecl32.03 


1998 


Sep 


27 


10:11 11:05 12:01 


2 


01:22:34 


10:18:22 


P 


ecl32.00 


1998 


Sep 


27 


10:17 11:11 12:07 





01:24:50 


08:54:49 


m 


ecl05.06 


1998 


Sep 


28 


2:21 3:37 4:37 


6 


22:48:57 


-01:02:56 


m 


ecl05.05 


1998 


Sep 


28 


2:27 3:43 4:43 


4 


22:51:14 


-02:26:12 


m 


ecl05.03 


1998 


Sep 


28 


2:33 3:49 4:49 


3 


22:53:32 


-03:49:22 


g 


ecl05.02 


1998 


Sep 


28 


2:46 3:55 4:55 


1 


22:55:50 


-05:12:33 


g 


ecl05.00 


1998 


Sep 


28 


2:52 4:01 5:01 





22:58:08 


-06:35:41 


g 


ecloe.oo 


1998 


Sep 


28 


2:58 4:07 5:07 





23:03:33 


-06:02:07 


g 


ecloe.oi 


1998 


Sep 


28 


3:04 4:13 5:13 


2 


23:01:13 


-04:39:05 


g 


ecl06.03 


1998 


Sep 


28 


3:10 4:19 5:19 


3 


22:58:54 


-03:15:59 


m 


ecl06.04 


1998 


Sep 


28 


3:16 4:25 5:25 


4 


22:56:36 


-01:52:54 


g 


ecl06.06 


1998 


Sep 


28 


3:22 4:31 5:31 


6 


22:54:18 


-00:29:46 


g 


e053m01.0 


1999 Feb 15 


10:00 10:54 11:50 


-1 


11:34:10 


01:47:35 


P 


e053_00.0 


1999 Feb 15 


10:07 11:00 11:56 





11:34:10 


02:47:35 


P 


e053p01.0 


1999 Feb 15 


10:12 11:06 12:02 


1 


11:34:10 


03:47:34 


P 
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Table 3 — Continued 



ID 


UT date 


UT times 










Objects 


e054p01.0 


1999 Feb 15 


10:18 11:12 12:08 


1 


11:37:49 


03:23:56 


P 




e054_00.0 


1999 Feb 15 


10:24 11:19 12:14 





11:37:49 


02:23:56 


P 




e054m01.0 


1999 Feb 15 


10:30 11:26 12:20 


-1 


11:37:50 


01:23:53 


P 




e055m01.0 


1999 Feb 15 


10:36 11:32 12:26 


-1 


11:41:29 


01:00:13 


P 




e055_00.0 


1999 Feb 15 


10:42 11:38 12:32 





11:41:29 


02:00:13 


P 




e055p01.0 


1999 Feb 15 


10:48 11:45 12:38 


1 


11:41:29 


03:00:13 


P 




e013p02.0 


1999 Feb 20 


3:14 4:33 5:45 


2 


09:00:36 


19:00:04 


m 




eOlSpOl.O 


1999 Feb 20 


3:20 4:40 5:51 


1 


09:00:36 


18:00:07 


m 




e013_00.0 


1999 Feb 20 


3:26 4:46 5:57 





09:00:36 


17:00:06 


m 




e013m02.0 


1999 Feb 20 


3:32 4:52 6:03 


-2 


09:00:36 


15:00:06 


g 




e016m01.0 


1999 Feb 20 


3:37 4:58 6:09 


-1 


09:12:45 


15:07:18 


m 




e016_00.0 


1999 Feb 20 


3:43 5:04 6:15 





09:12:45 


16:07:17 


g 




eOlGpOl.O 


1999 Feb 20 


3:49 5:10 6:21 


1 


09:12:45 


17:07:18 


m 




e016p02.0 


1999 Feb 20 


3:55 5:16 6:27 


2 


09:12:45 


18:07:18 


m 




e023p02.0 


1999 Feb 20 


4:01 5:22 6:33 


2 


09:40:39 


15:54:37 


m 




e023p01.0 


1999 Feb 20 


4:07 5:28 6:39 


1 


09:40:39 


14:54:36 


g 




e023_00.0 


1999 Feb 20 


4:13 5:33 6:45 





09:40:39 


13:54:37 


m 




e023m01.0 


1999 Feb 20 


4:19 5:39 6:51 


-1 


09:40:39 


12:54:35 


g 




e032m01.0 


1999 Feb 20 


6:57 8:10 9:21 


-1 


10:15:38 


09:47:45 


g 


1999 DEg 


e032_00.0 


1999 Feb 20 


7:03 8:16 9:27 





10:15:38 


10:47:45 


g 




e032p01.0 


1999 Feb 20 


7:09 8:22 9:33 


1 


10:15:38 


11:47:45 


g 


1999 DFg 


e032p02.0 


1999 Feb 20 


7:15 8:28 9:39 


2 


10:15:38 


12:47:44 


g 




e040m01.0 


1999 Feb 20 


7:46 8:58 10:22 


-1 


10:45:58 


06:50:10 


g 




e040_00.0 


1999 Feb 20 


7:52 9:03 10:28 





10:45:58 


07:50:11 


g 




e040p01.0 


1999 Feb 20 


7:58 9:09 10:34 


1 


10:45:58 


08:50:11 


g 




e040p02.0 


1999 Feb 20 


8:04 9:15 10:40 


2 


10:45:58 


09:50:11 


g 




e047m01.0 


1999 Feb 20 


10:46 11:36 12:12 


-1 


11:12:03 


04:08:46 


g 




e047_00.0 


1999 Feb 20 


10:52 11:42 12:18 





11:12:03 


05:08:45 


m 




e047p01.0 


1999 Feb 20 


10:58 11:48 12:25 


1 


11:12:03 


06:08:46 


g 




el002p05 


2000 Mar 03 


3:06 4:22 5:31 


5 


09:24:33.7 


20:18:02.1 


g 




el002p04 


2000 Mar 03 


3:19 4:28 5:37 


4 


09:24:33.7 


19:17:02.0 


g 




el002p03 


2000 Mar 03 


3:25 4:35 5:44 


3 


09:24:33.7 


18:16:01.9 


g 
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Table 3 — Continued 



ID 


UT date 


UT times 


13 * 


a 


5" 




Objects 


el002p02 


2000 Mar 03 


3:31 4:41 5:50 


2 


09:24:33.7 


17:15:01.8 


g 




el002p01 


2000 Mar 03 


3:37 4:47 5:57 


1 


09:24:33.7 


16:14:01.6 


g 




el002_00 


2000 Mar 03 


3:44 4:53 6:03 





09:24:33.7 


15:13:01.5 


g 




el002m01 


2000 Mar 03 


3:50 5:00 6:10 


-1 


09:24:33.7 


14:12:01.4 


g 




el002m02 


2000 Mar 03 


3:57 5:06 6:15 


-2 


09:24:33.7 


13:11:01.3 


g 




el002m03 


2000 Mar 03 


4:03 5:12 6:22 


-3 


09:24:33.7 


12:10:01.2 


g 


2000 EEi73 


el002m04 


2000 Mar 03 


4:09 5:18 6:28 


-4 


09:24:33.7 


11:09:01.0 


g 




el002m05 


2000 Mar 03 


4:15 5:25 6:34 


-5 


09:24:33.7 


10:08:00.9 


g 




el020m05 


2000 Mar 03 


6:42 7:51 9:00 


-5 


10:40:33.6 


03:17:43.3 


g 




el020m04 


2000 Mar 03 


6:48 7:57 9:06 


-4 


10:40:33.6 


04:18:43.4 


g 




el020m03 


2000 Mar 03 


6:54 8:03 9:13 


-3 


10:40:33.6 


05:19:43.5 


g 




el020m02 


2000 Mar 03 


7:00 8:10 9:19 


-2 


10:40:33.6 


06:20:43.6 


g 




el020m01 


2000 Mar 03 


7:07 8:16 9:26 


-1 


10:40:33.6 


07:21:43.8 


g 




el020_00 


2000 Mar 03 


7:13 8:22 9:32 





10:40:33.6 


08:22:43.9 


g 




el020p01 


2000 Mar 03 


7:19 8:28 9:38 


1 


10:40:33.6 


09:23:44.0 


g 




el020p02 


2000 Mar 03 


7:24 8:35 9:45 


2 


10:40:33.6 


10:24:44.1 


g 




el020p03 


2000 Mar 03 


7:32 8:41 9:51 


3 


10:40:33.6 


11:25:44.2 


g 




el020p04 


2000 Mar 03 


7:38 8:47 9:57 


4 


10:40:33.6 


12:26:44.4 


g 




el020p05 


2000 Mar 03 


7:45 8:54 10:03 


5 


10:40:33.6 


13:27:44.5 


g 




el041m05 


2000 Mar 03 


10:11 10:56 11:55 


-5 


12:05:04.3 


-05:37:58.9 


m 




el041m04 


2000 Mar 03 


10:17 11:03 12:01 


-4 


12:05:04.3 


-04:36:58.8 


g 




el041m03 


2000 Mar 03 


10:23 11:09 12:08 


-3 


12:05:04.3 


-03:35:58.6 


g 




el003m05 


2000 Mar 04 


2:28 3:44 4:47 


-5 


09:28:54.5 


09:47:19.9 


m 




el003m04 


2000 Mar 04 


2:34 3:50 4:54 


-4 


09:28:54.5 


10:48:20.0 


P 




el003m03 


2000 Mar 04 


2:40 3:57 5:06 


-3 


09:28:54.5 


11:49:20.1 


m 




el003m01 


2000 Mar 04 


2:59 4:03 5:13 


-1 


09:28:54.5 


13:51:20.4 


m 




el003_00 


2000 Mar 04 


3:05 4:10 5:20 





09:28:54.5 


14:52:20.5 


m 




el003p01 


2000 Mar 04 


3:13 4:16 5:26 


1 


09:28:54.5 


15:53:20.6 


m 




el003p02 


2000 Mar 04 


3:19 4:22 5:32 


2 


09:28:54.5 


16:54:20.7 


m 




el003p03 


2000 Mar 04 


3:25 4:29 5:38 


3 


09:28:54.5 


17:55:20.8 


m 




el003p04 


2000 Mar 04 


3:32 4:35 5:45 


4 


09:28:54.5 


18:56:21.0 


m 




el003p05 


2000 Mar 04 


3:38 4:41 5:51 


5 


09:28:54.5 


19:57:21.1 


g 
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Table 3 — Continued 



ID UT date UT times /3 a ^ (5 ^ 9'^ Objects 



el021m04 


2000 Mar 04 


5:58 6:55 i 


3:14 


-4 


10:44:39.7 


03:54:03.5 


m 


el021m03 


2000 Mar 04 


6:04 7:01 i 


3:20 


-3 


10:44:39.7 


04:55:03.6 


m 


el021m02 


2000 Mar 04 


6:10 7:07 i 


3:27 


-2 


10:44:39.7 


05:56:03.7 


m 


el021m01 


2000 Mar 04 


6:16 7:25 i 


3:33 


-1 


10:44:39.7 


06:57:03.8 


m 


el021_00 


2000 Mar 04 


6:23 7:32 i 


3:39 





10:44:39.7 


07:58:04.0 


m 


el021p01 


2000 Mar 04 


6:29 7:40 i 


3:45 


1 


10:44:39.7 


08:59:04.1 


P 


el021p02 


2000 Mar 04 


6:35 7:46 i 


3:52 


2 


10:44:39.7 


10:00:04.2 


P 


el021p03 


2000 Mar 04 


6:41 7:52 i 


3:58 


3 


10:44:39.7 


11:01:04.3 


P 


el035p05 


2000 Mar 04 


9:26 10:17 11:10 


5 


11:41:11.8 


07:07:06.2 


m 


el036p05 


2000 Mar 04 


9:32 10:23 11:17 


5 


11:45:11.1 


06:41:15.9 


m 


el037p05 


2000 Mar 04 


9:39 10:29 11:23 


5 


11:49:10.1 


06:15:24.8 


P 



Note. — This table lists fields imaged with the KPNO 8k Mosaic camera. Fields were 
imaged in triplets, with UT times given for each image. KBOs found are listed after the field 
of discovery. The complete version of this table is in the electronic edition of the Journal. The 
printed edition contains only a sample. 

^J2000 ecliptic latitude, degrees 
^J2000 right ascension, hours 
'^J2000 declination, degrees 

■^Seeing category: g, m, p represent the good (< 1.5 arc sec), medium (> 1.5 arc sec and 
< 2.0 arc sec), and poor (> 2.0 arc sec) seeing cases, respectively. The efficiency functions for 
each of these cases are presented in Table § 
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Table 4. KPNO Survey Efficiency 





Good 


Medium 


Poor 


Global 


Median PSF FWHM ["] 


1.3 


1.6 


2.1 


1.5 


PSF FWHM Range ["] 


1.0-1.4 


1.5-1.9 


2.0-2.4 


1.0-2.4 


6max 


0.88 


0.88 


0.80 


0.85 




21.4 


21.2 


20.8 


21.1 


a 


0.5 


0.5 


0.4 


0.6 


Fields Imaged 


57 


61 


53 


171 


Total Area [deg^] 


55 


58 


51 


164 


KBOs 


3 








3 


Centaurs 


1 








1 



Note. — Survey efficiency for three seeing cases and "Global" 
efficiency, representing the sky-area weighted mean of the three 
efficiencies, designed to represent the properties of the KPNO data 
set in its entirety. 
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Table 5. The KPNO Objects 





1998 


1999 DEg 


1999 DFg 


2000 EE173 


Discovery Conditions 










Date 


1998 Sep 26 


1999 Feb 20 


1999 Feb 20 


2000 Mar 03 


R [AU] 


38.248 


32.549 


39.765 


23.593 


A [AU] 


37.264 


33.537 


38.776 


22.678 


a' [deg] 


0.3 


0.0 


0.0 


0.9 


rriYi 


21.2 ±0.1 


19.7 ±0.1 


21.6 ±0.1 


22.2 ±0.1 


mR(l,l,0) 


5.4 


4.5 


5.7 


8.4 


Diameter [km] 


480 


710 


420 


120 


a [AU] 


37.900 ± 0.003 


56.4 ± 0.5 


46.400 ± 0.053 


49.02 ± 0.06 


e 


0.0437 zb 0.0004 


0.435 ± 0.009 


0.144 ± 0.004 


0.5387 ± 0.0007 


i [deg] 


4.605 ± 0.0005 


7.623 ± 0.002 


9.819 ±0.005 


5.950 ± 0.0004 


(node) [deg] 


192.097 ± 0.002 


322.983 ± 0.003 


334.856 ± 0.002 


294.070 ± 0.002 


uj (peri) [deg] 


264.156 ± 0.145 


155.7 ±1.7 


169 ± 11 


235.05 ± 0.13 


T (time of peri) [JD] 


2474593 ± 36 


2445972 ± 222 


2449415 ± 2573 


2454160 ± 10 


epoch 


2451071.7 


2451229.8 


2451229.8 


2451606.7 


perihelion 


36.33554 


31.83661 


39.70701 


22.62496 



Note. — Discovery conditions of the objects found in this survey. R is the heliocentric distance. 
A is geocentric distance, a' is phase angle. mR is the red magnitude at discovery. mR( 1,1,0) 
is red magnitude at zero phase angle, with heliocentric and geocentric distance = 1 AU. The 
diameter is derived assuming a geometric albedo of 4%. Orbital elements and associated errors 
were computed using the procedure of Bernstein and Khushalani (2000) with data provided by 
the Minor Planet Center. 

^This object was discovered by Spacewatch on Sep 23, 1998 and serendipitously detected in 
our data. 
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Table 6. Classical KBO Model Assumptions 



Quantity 


Value 


Distribution 


Description 


a [AU] 


41< a < 47 


n{a)da ~ a}~^da 


semimajor axis distribution 


pa. 


2 




semimajor axis power 




q' < 37.0 




perihelion criterion 


i 




Gaussian 


inclination distribution 


«l/2ldGg] 


20.0 




inclination half-width 


r [km] 


25 < r < 16000 


n{r)dr ~ r~'^dr 


radius distribution 


q 


fitted 




radius distribution power 


iVcKBOs(^ > 100 km) 


fitted 




number of CKBOs with D > 100 km 


A^bins 


50 


logarithmic 


number of radius bins 



^For circular orbits, p represents the decrease in ecliptic plane surface density Sed as a function of 
heliocentric distance R, Heel ~ R~^- 

^Eccentricity assumed to be uniform between and 1 consistent with perihelion criterion. 
'^Best-fit values from Trujillo, Jewitt Sz Luu (2001). 
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Table 7. Objects Brighter than 20 R 



1996 GQ21 1999 DE9 2000 EB173 



Discovery Conditions 



Date 


1996 Apr 12 


1999 Feb 20 


2000 Mar 10 


R [AU] 


38.523 


32.549 


29.940 


A [AU] 


37.521 


33.537 


29.036 


a' [deg] 


0.1 


0.0 


0.8 


m 


20.5 V 


19.7 R 


19.3 R 


m(l,l,0) 


4.7 V 


4.5 R 


4.5 R 


Diameter [km] 


720 


710 


780 


a [AU] 


94.06 ±0.10 


56.4 ± 0.5 


39.22 ± 0.03 


e 


0.5955 ± 0.0007 


0.435 ± 0.009 


0.2697 ±0.0011 


i [deg] 


13.357 ± 0.0003 


7.623 ± 0.002 


15.476 ± 0.0004 


(node) [deg] 


194.184 ± 0.0003 


322.983 ± 0.003 


169.331 ± 0.001 


uj (peri) [deg] 


353.992 ±0.3 


155.7 ± 1.7 


66.8 ± 0.2 


(time of peri) [JD] 


2447365 ± 54 


2445972 ± 222 


2456944 ± 22 


epoch [JD] 


2450185.8 


2451229.8 


2450182.7 


source 


Larsen et ah (2001) 


This Work 


Ferrin et al. (2001) 



Note. — Discovery conditions of ah KBOs brighter than 20 R discovered in pubhshed 
surveys (Pluto excepted). R is the hehocentric distance. A is geocentric distance, m is 
the magnitude at discovery. m(l, 1, 0) is magnitude at zero phase angle, with geocentric 
and heliocentric distance = 1 AU. The diameter is derived assuming a geometric albedo 
of 4%. Orbital elements and associated errors were computed using the procedure of 
Bernstein and Khushalani (2000) with data provided by the Minor Planet Center. 



